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ABSTRACT 

In  this  paper  we  inv('stigate  the  nonlinear  development  of  the  most  nnstable 
Gdrth'r  mod('  within  a  general  three'  diiiK'nsional  honndary  layer  upon  a  suital)ly 
concave  surface.  The  structure  of  this  mode  was  first  identified  by  Denier,  Hall 
Seddougui  (1991)  who  demonstrated  that  the  growth  rate  of  this  instability  is  0(G^  ) 
wlu'ie'  G  is  tlu'  Gdrth'r  numbe'r  (taken  to  Ix'  large'  here),  which  is  effectively  a  measure' 
e)f  the'  e-urvature'  e)f  the'  surface'.  Previems  re'.seare'hes  have  elescribe'd  the  fate'  e)f  the 
most  unstable'  me)ele'  within  a  twe)  elinu'iisieMial  boundary  laye'r.  De'iiier  Hall  (1992) 
elise-usse'el  the'  fully  nonliiK'ar  ele've'le)pme'nt  e)f  the'  vente'x  in  this  e-ase  and  shenveel  that 
the'  ne)idine’arity  e-ause's  a  bre'akele)wn  of  the  fle)w  structure. 

The’  e'ffe'ct  e)f  e-re).ssfle)W  anel  unste’aeline'ss  upem  an  infinitesimal  unstable  me)ele'  was 
e’lue-ielate’el  by  Bassenn  ^  Hall  (1991).  The'y  ele'me)nstrate'el  that  crossflow  te'iids  te) 
stabilise'  the-  me).st  unstable’  Gdrth'r  nmele,  anel  fe)r  e-e’rtain  crossflenv/freepie'iicy  combi- 
natienis  the'  Ge'irth'r  nie)eh’  may  be  maeh'  ne’utrally  stable.  The'.se  veute'X  e'eniflgurations 
naturally  h'liel  the'inse'lves  to  a  we’akly  neudine'ar  stability  analysis;  we)rk  whiedi  is  ch’- 
scribe'el  in  a  i)re’vie)us  artie-h’  by  the  pre'se'nt  authens.  He're  we  e’xte’iid  the  ieh'as  of 
De'iiie'r  anel  Hall  (1992)  te)  the’  thre’e'  elime’nsie)nal  beuuidaiy  layer  problem.  It  is  fe)unel 
that  the’  nume-rical  se)lutiem  e)f  the  fedly  ne)nline’ar  eepiations  is  be’st  ce)nelucteel  using  a 
nie’the)el  whie'h  is  e’sse’iitiaily  an  aelaptie)n  e)f  that  utilise’el  by  De’iiie’r  anel  Hall  ( 1992).  The 
influe'iice’  e)f  cre)ssfle)W  anel  unste’aeliiu’ss  upein  the  breakehiwn  of  the  fleiw  is  elescribeel. 
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§1  Introduction 


The  aim  of  this  article  is  to  further  our  \inderstancliugof  the  effects  of  unsteadiness 
and  crossflow  ui>on  the  fully  nonlinear  development  of  unstable  Gdrtl(>r  modes.  The 
initial  derivation  of  the  governing  ecpiatioiis  for  these  modes  was  given  by  Gortler 
(1940)  whos('  results  were  modified  by  Harnmerlii:  (1956).  These  early  works  ignored 
non})arallel  ('ff('cts  pn'sent  due  to  boundary  layer  growth  and  Smith  (1955)  added 
t('rms  in  an  attempt  to  rectify  this  deficiency.  Until  recent  years  it  was  unclear  as  to 
the  importance  of  the  nonparallel  terms;  this  question  was  resolved  by  the  results  of 
Hall  (1982  a,b,  1983).  In  this  series  of  papers  Hall  showed,  using  a  combination  of 
asymptotic  and  numerical  technicpies,  that  for  order  one  wavenumber  vortices  there  is 
no  unicpie  m'utral  linear  stability  curve.  More  precisely,  the  stability  characteristics  of 
such  wav(’numb('r  modes  are  entirely  dependent  u])on  the  initial  form  and  location  of 
the  disturbance.  However  for  small  wavelength  vortices  a  unicpie  neutral  ciirvc'  does 
exist  and  on  this  curve  the  vort('x  wavenumber  k  is  0{G^)  where  G  is  the  (assumed 
larg(')  Giirth'r  number.  For  further  reading  concerning  the  develo]unent  of  the  stability 
theory  for  Gfirtler  vortices  the  reader  is  referred  to  Hall  (1990). 

The  most  unstable  Gortler  mode  (i.e.  that  infinitesimally  sized  vortex  which  has 
the  largest  growth  rate)  was  ol)tained  by  Denier,  Hall  &  Seddougui  (1991)  and  Tim- 
oshin  (1990).  By  considering  the  stability  properties  of  the  essentially  viscous  modes 
of  wavenumber  0(G^ )  together  with  those  of  the  inviscid  modes  of  0(1)  wavelengths 
it  was  possilih'  to  identify  an  int('rin('d’atc  wavenumber  ix'gime  in  which  the  vortc'x 
growth  rat('  is  larg(’st.  These  unstable  modes  are  confined  to  a  region  of  dejitli  of 
()(G~^ )  and  possess  growth  rates  of  size  0(G^ ).  The  stability  properties  of  0(0^ ) 
wavenumlier  vortices  are  deduced  by  solving  a  sixth-order  ordinary  differential  system 
with  appropriate  boundary  conditions  and  the  solution  of  this  systcnn  reveals  that  the 
unicjue  most  unstable  mode  has  wavenumber  k  =  0.476G^  and  growth  rate  0.312G's. 

To  date  tlu're  has  bei'ii  relatively  little  attention  paid  to  the  nonlinear  staliility 
jiroperties  of  Gortler  vortices.  Perhaps  the  first  work  was  jierfonned  by  Aihara  (1976) 
who  attempted  to  descrilie  the  nonlim'ar  evolution  of  Gortler  vortices  using  paralh'l 
argunu'nts.  Later  calculations  by  Hall  (1988)  remedied  these  defects  and  showed  that 
as  ()( 1 )  wav<'number  modes  evolve  downstream  so  the  energy  of  tlu'  flow  concentrati's 
itself  in  the  fundamental  and  mean  flow  correction.  This  sugg(’sts  far  downstream  tht' 
flow  can  l)e  adecjuately  d(’scribed  l)y  a  mean  fi('ld/first  harmonic  structure  and  such  a 
configuration  was  elucidated  for  short  wavelength  modes  using  Ijoth  wt-akly  nonlinear 
and  fully  nonlinear  approaclu’s  l)y  Hall  (1982  !>)  and  Hall  Lakin  (1988)  resp('ctiv('ly. 

A  fuller  description  of  the  nonlinear  stability  of  Gortler  iikkIcs  may  be  found  in 
Denier  &  Hall  (1992).  In  that  paper  the  authors  argm'd  that  in  a  numbt'r  of  practical 
syst('ms,  ('sp('cially  wlu're  significant  curvature  occurs  such  as  tlu’  cas('  of  flow  ov<'r 
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turbine  Ijlades,  one  would  expect  the  small  localized  surface  imperfections  may  well 
trigger  the  most  unstable  linear  Gdrtler  mode.  (This  conclusion  relies  on  the  result 
of  the  linear  receptivity  theory  given  by  Denier,  Hall  &:  Seddougui  (1991).)  This 
motivated  a  careful  study  of  nonlinear  evolution  of  the  most  unstable  mode  which  was 
tackled  in  the  following  way.  Denier  &  Hall  (1992)  took  an  arbitrary  form  of  initial 
disturbance  at  a  typical  streamwise  location,  say  .r  =  1.  By  integrating  the  linear 
form  of  the  vortex  equations  over  a  long  distance,  up  to  j-  =  101,  a  flow  profile  was 
obtained  which  is  dominated  by  the  most  unstable  mode  described  by  Denier  et  al 
(1991).  From  x  —  101  onwards  the  full  nonlinear  equations  replaced  the  linear  ones 
and  the  most  unstable  vortex  mode  was  marched  further.  Typically  8  or  16  harmonics 
of  the  fundamental  were  retained  during  this  calcidation.  It  was  found  that  at  a  critical 
point  the  flow  contains  a  regioti  of  reverse  flow  and  the  analysis  is  then  no  longer  valid. 
Denier  &  Hall  (1992)  interpreted  this  breakdown  as  being  responsible  for  the  vortices 
moving  away  from  the  wall  and  into  the  core  of  the  boundary  layer. 

The  effect  of  crossflow  and  unsteadiness  on  the  most  unstable  Gdrtler  mode  was 
discussed  by  Bassom  &  Hall  (1991).  The  primary  result  arising  from  this  work  was 
the  demonstration  that  a  relatively  small  crossflow  could  completely  stabilise  the  most 
unstable  mode.  Additionally,  by  allowing  for  vortex  unsteadiness,  it  was  shown  that 
suitable  combinations  of  crossflow,  vortex  frequency  and  wavenumber  could  lead  to 
neutrally  stable  configurations.  A  weakly  noiilinear  stability  analysis  pertaining  to 
such  configurations  was  conducted  by  Bassom  &  Otto  (1992)  who  derived  classical 
‘Stuart- Watson’  (1960)  type  evolution  equations  for  near-neutral  inodes.  They  con¬ 
cluded  that  the  weak  nonlinearity  has  a  stabilising  effect  and  derived  equations  for 
the  supercritical  equilibrium  amplitudes. 

The  results  of  Denier  &  Hall  (1992)  and  Bassom  &:  Hall  (1991)  provide  the  mo¬ 
tivation  for  the  current  study.  Within  a  two-dimensional  boundary  layer  the  effect  of 
nonlinearity  on  the  most  unstable  mode  tends  to  lead  to  a  finite-distance  breakdown 
whereas  crossflow  appears  to  stabili.se  the  flow.  With  these  two  mechanisms  tending 
to  have  opposite’  e’ffects  it  is  clear  that  in  many  practical  situations,  where  tlirc'e 
fliniensionality  is  undoulite’dly  important,  it  is  of  great  interest  to  determine  which 
of  these  two  conflicting  behaviours  dominate.  We  attempt  to  answer  this  question  by 
considering  the  full  nonlinear  vortex  equations  and  employ  numerical  techniques  which 
are  similar  to  those  used  in  Denier  &  Hall  (1992)  but  modified  in  certain  ways  (detaile<l 
lat('r).  These  im])rovements  significantly  speed  up  the  computations  and  allow  us  to 
obtain  a  greater  range  of  results  than  those  found  by  Denier  Hall  (1992). 

Tlu’  structure  of  the  r<’mainder  of  this  article  is  as  follows;  in  section  2  th('  fully 
nonlinear  equation'^  are  deriv(’d  and  a  brief  descrii)tion  of  the  munerical  procedures 
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used  are  outlined  in  section  3.  Details  of  the  results  are  presented  in  section  4  and  in 
section  5  some  brief  conclusions  are  drawn. 


§2  Formulation 

As  in  Hall  (1985)  we  consider  a  boundary  layer  flowing  over  the  cylinder  y  =  0, 
—oo  <  5  <  oo  where  the  f-axis  is  a  generator  of  the  cylinder,  y  measures  the  distance 
normal  to  the  surface  and  x  denotes  the  distance  along  the  surface.  The  Reynolds 
number  Be  ^md  Gdrtler  number  G  are  defined  by 


G  =  2rU, 


where  Uq  is  a  typical  flow  velocity  in  the  i-direction,  L  is  a  characteristic  streamwise 
lengthscale  and  u  is  the  kinematic  viscosity  of  the  fluid.  Furthermore  the  curvature  of 
the  cylinder  is  supposed  to  be  (x)  and  with  the.se  definitions  6  =  T/ft,  where  h  is 
a  typical  radius  of  curvature  of  the  cylinder  (Gortler  vortices  are  typically  observed  in 
flows  over  concave  surfaces  which  corresponds  to  the  choice  xo  >  0  ).  The  Reynolds 
number  is  supposed  to  lie  large  whilst  6  is  sufficiently  small  so  that  as  (ii  — >  0  the 
parameter  G  is  fixed  and  is  of  order  one.  The  basic  three-dimensional  boundary  layer 
is  taken  to  be  of  the  form 


u  =  17o  (n  ( -Y,  r ), Re  ( X, F ), Re  ^ A* «.  (X,  F ))( 1  +  O  (i?r ^ ) ) 

i 

where  X  =  x/ L  and  F  =  yRj  j L  and  the  cro.ssflow  parameter  A*  is  of  order  one. 

i 

It  is  convenient  to  define  the  scaled  spanwise  coordinate  Z  =  zRf  / L  and  let  T  be 
the  temporal  varialjle  scaled  on  L/Uq.  The  basic  velocity  profile  is  perturbed  by  the 
quantity 

(f^(r,X,F,Z),Re"V(T,X,F,Z),Re"%(r,X,F,Z))  , 

_  I 

and  the  pressure  field  by  Re  ^R(r,  X,  F,  Z).  Substituting  this  flow  form  into  the 
continuity  and  Navier  Stokes  equations  yields  the  system 


Ux  +  Vy  +Wz  =  0, 


(2.1n) 


-Ut  +  Uyy  +Uzz  -  uyV  -  uUx-TixU  -  vUy  -  A*Tcf7, 

=  UUx  +  VUy  +  WUz, 

—  Vp  +  Vyv  +  Vzz  ~  G\uU  —  Py  —  ‘flVx  —  f’xU  —  vVy  —  cy  V  —  \*U'V~ 


(2.16) 


G  (2.1c) 

UVx  +  FFy  +  WVz  +  -\U\ 
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(2.  If/) 


—  Wt  +  Wyy  +  Wzz  —  Pz  —  uW X  —  \*tvxU  —  vWy  —  A*V"?/)v  —  \*wW, 

=  UWx  +  VWy  +  WWz, 

_  1 

where  terms  of  relative  order  0(Re  ^ )  have  been  neglected.  It  is  worth  noting  at  this 
point  that  the  linearised  system  studied  by  Bassoni  k  Hall  (1991)  is  oljtained  liy  setting 
the  right-hand  sides  of  (2.1  b-d)  to  zero  whereas  the  nonlinear  ecpiations  examined 
by  Denier  k  Hall  (1992)  to  determine  the  development  of  steady  nonlinear  vorticf's  in 
two  dimensional  boundary  layers  can  be  retrieved  I>y  setting  A*  and  dr  equal  to  zcuo. 

We  now  invoke  the  scalings  proposed  by  Denier  et.  al.  (1991)  who  flemonstrated 
that  in  high  Gdrtler  number  flows  the  most  unstable  vortices  have  O(G^)  wavenuml)ers 
and  are  confined  to  a  layer  of  thickness  of  0{G~^)  adjacent  to  the  cylinder.  These 
modes  have  a  spatial  growth  rate  of  0{G^)  and  we  use  the  residts  of  Bassom  k  Hall 
(1991)  who  illustrated  that  the  three  dimensionality  of  the  l^asic  flow  significantly 
affects  the  two  dimensional  staliility  results  once  the  scaled  crossflow  parameter  A* 

3  A 

becomes  0(G^  ).  Therefore  it  is  convenient  to  define  the  0(1)  crossflow  parameter  A 

by 

A*  =  Gti  (2.2o) 

To  reflect  the  fact  that  the  vortices  are  confined  to  a  region  of  depth  0{G~^  )  adjacent 
to  the  cylind('r  we  introduce  the  0(1)  boundary  layer  coordinate  y  defined  by 


y  =  on; 


{2.2b) 


and  in  this  layer  the  basic  flow  may  l)e  expanded  as  a  Maclaurin  series  of  the  form 


u  =  0-  =  //n(A')y  +  iG-t/fr2(A)?/  +  j^G-^iniX),/  +  ■■■ 
w  =  G~^fi-2\{X)y  +  ^G~^/i22iX)i/  +  ^G~^-  fi2:i(X)y^  H - 


(2.2r) 

(2.2f/) 


To  determine  the  form  of  the  vortf'x  disttirbance  we  appeal  to  the  findings  of  Bassom 
A  Otto  (1992)  who  identifii'd  the  crucial  perturl)ation  size  at  which  the  governing 
cfjuations  l)ecome  fully  nonlinear  (although  these  authors  inach'  no  attempt  to  solv<' 
these  fully  nonlinear  forms).  The  disturbance'  forms  are  then 


G  =  G-^  (u^  +  G-^u^  +G-h^2+  -  ).  V  =  G^  (v), +  G-h;  +G-h; +  •••). 

{2.3a. h) 

ir  =  G ^  (h  0  +  G-  s  ir,  +  G- ^  H',  +  •  •  •)  .  P  =  G s  (Po  +  G-  ^  P,  +  G- ^  P,  +  ■  ■  ■)  . 

(2.3r.d) 

where  fO,  \’().  H  0.  Pi.  Pi .  ■  .  .  etc.  are  all  functions  of  AT  y.  tlu'  temporal  and  the  sjiau- 
wise  varial)l('.  It  is  now  convenient  to  imph'me'nt  the  residt  of  Bassom  k  Hall  (1991) 
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that  the  leading  order  liehaviour  in  the  downstream  coordinate  is  purely  oscillatory 
and  we  do  this  by  introducing  the  coordinate  and  temporal  variations  given  Ijy 


=  Z  =  +  T=G-^t, 


(2.4n) 


where  /i  =  A/<2i//iii-  This  then  implies  that  the  streamwise  and  spanwise  derivatives 
become 

d  d  1  d  d  d  ,  , 

The  desired  governing  equations  are  obtained  by  substituting  (2.2)-(2.4)  into  equations 
(2.1).  Leading  order  terms  in  the  momentum  equations  yield  that  I4^o  =  f^Uo-  Next 
order  terms  in  equation  (2.1  b)  give  the  first  relationship 


/  y'^a  d  d  d 

~WT~ 


Uo  —  i^t^x  1  Vo  =  RHSi , 


(2.5n) 


where  the  precise  form  of  RHSi  is  given  presently  and  cv  =  A/i22  —  l^yzx-  A  second 
equation  is  derived  by  following  a  procedure  similar  to  that  described  Ijy  Bassom  Hall 
(1991).  By  considering  the  second  order  terms  in  the  y  and  c  momentum  equations 
(2.3),  eliminating  the  pressure  by  cross-differentiation  and  applying  the  continuity 
equation  it  is  a  routine  but  lengthy  task  to  obtain 


dVo 

dz 


dHro 

dz^ 


=  RHS2 


^  dz^  2  dz  dT  )  Vay2  dz^ ) 

(2.5/;) 

where  again  we  shall  specify  RHS2  shortly.  The  forms  of  RHS2  includes  reference  to 
the  combination  Wi  —  PU\  which  therefore  needs  to  be  expressed  in  terms  of  qiiantities 
with  subscript  zero.  This  is  best  accomplished  by  integrating  the  continuity  equation 
to  give 

\dUo  ,  aVo' 


-du,  =<^o(x 


,y)-j 


+ 


dz. 


(2.5c) 


dx  dy 

The  aim  of  this  article  is  to  consider  the  nonlinear  evolution  of  modes  which  are  periodic 
in  the  spanwise  direction,  with  fundamental  wavenumber  (xo/tn)  ^  k  say,  and  periodic 
in  time.  Now  it  is  advantageous  to  introduce  scalings  first  proposed  by  Bassom  <Sc  Hall 
(1991).  It  is  found  that  if  we  transform  according  to 


dy  (xo/Zf,)"  kdy,  d,  (xo//n)'  kd,.  dr  ^  X^hn^x, 


Uq  =  Xo 


=  Xo/'n'/» 


=  \  0  r  1 V,  d-r  X  0  /'  n  ^7' ' 

A  ifli  1/122  —  I'ZX/lxz) 


f2.6a  —  c) 
{2.6d-f) 
{2.6y-h) 


thoii  th<'  transfornu'd  ('quations  (2.5  a,b)  are  independent  of  the  i^articnlar  boundary 
layer  innh'r  eon.sideration.  Therefore  the  ensiling  results  are  jiotentially  relevant  to 
a  wid('  class  of  liasic  configurations.  In  this  nondiniensionalised  coordinate  system 
w’v  restrict  oursidves  to  flow  quantities  the  fundamental  mode  having  sjianwisi'  and 
ti'inporal  variation  given  liy  exj)  [i  (~  +  ^0]’  represents  the  nondimensional 

frt'quency  of  the  fundamental  vortex  component. 

Ui>on  making  the  transformations  (2.6)  the  leading  order  vortex  equations  (2.5 
a.b)  become 


r  .  O' 

Ol)^  ^  Oz^ 


,  2A  dV 


y 


F  dzOy  ~  P  dxdy 


P  dz  P  dz^  P  dz^ 


.  (2.7r;) 


wIk’H'  till'  operator  L  is  defiiu'd  according  to 


(2.7/1) 


\ij^  dir  1  dii'  y  dij' 

^  ^  "  Pdf  ~  Pdf' 


(2.7c) 


and 


s.n  =  r^,,r^  +  4 

dj-  dy  Dz 


1  f\dU 

A-/[ar 


+  k 


Oy . 


dz 


+  v 


-2 


=r 


d 


d.r 

+  V'A 


d  ■ 

d 


dy 

-ll 


II 
HI 


\dU  dV 
dx  dy  J 


du  ,  dV 

+  k-~ 


dz 


ofL  2  ^ 

dx  ""  dy  J 


dx  dy  \  J 
dV 


dz 


dV\ 

0, 


^.dU  ,^^dU  dU 


o 


-\L 


dU  dV 

dx  '^^'dy 


dz 


{2.7d) 


(2.7c) 


(2.7/) 


To  comjilete  tlu'  governing  equations  it  is  lu'cessary  to  determine  d.  This  is  achii'ved 
by  noticing  that  U'l  —  dU\  satisfies 
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Identifying  tlu'  component  of  this  <'quatiou  imlejxuident  of  c  yields 
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To  close  the  system  requires  sjH'cificatiou  of  appropriate  l)oun(lary  coii(litioiis.  Clearly 
we  require  the  velocity  coinponeiits  U.  V'  to  vanish  ou  ,(/  =  0  (and  In'  continuity  so 
must  dVjdy).  Additionally  tlu'  mean  flow  terms  <3  =  0  on  y  =  0  and  in  order  that 
the  disturhance  he  confined  to  the  boundary  layer  we  fUmiand  that  the  str('amwis(' 
velocity  U  tends  to  some  function  of  r  as  y  oc. 

We  notice  that  (’quations  (2.7)  are  the  appropriate  geiK'ralisations  of  those  solved 
by  Denier  and  Hall  for  the  most  unstable  nonlinear  vortex  within  a  two  dimensional 
boTuidary  layer  (tludr  equations  are  r<'cover<'d  l>y  setting  A  =  <;5  =  0  and  setting  d'j  =0). 
This  allowi'd  us  t(i  comi)are  numerical  results  against  their  i^reviously  published  ont's 
as  a  check  of  our  numerical  methods. 

§3  Numerical  Methods 

The  methods  employed  to  solve  system  (2.7)  are  similar  in  ('ssence  to  those'  use'd 
in  Denier  <k:  Hall  (1992).  However  it  was  found  necessary  to  introduce  a  number  of 
amendnu'nts  to  their  code  in  order  to  speed  up  the  computations  which  in  turn  had 
the  benefit  of  allowing  an  increa.sed  number  of  results  to  be  obtaiiu'd. 

Denier  Hall  (1992)  solve'd  the*  two  (limensional  counte'iparts  of  ('(piations  (2.7) 
as  follov's.  They  decompo.sed  each  of  the  flow  velocities  C.V'.  IT  into  tlu'ir  Fouric'r 
comi)onents  and  rewrote  the  governing  equations  in  terms  of  these  components.  By 
utilising  a  scheme  based  upon  that  implemented  by  Hall  (19SS)  they  obtaiiK'd  finite- 
difFer<'nc('  ('quations  for  the  mean  flow  and  harmoni<  t('rins.  These'  ('(luations  contain 
only  one  streamwise  derivative  and  a  straightforward  method,  based  upon  solving  one 
tridiagonal  and  one  pentadiagonal  syste'in.  may  be  us(’d  to  march  the'  .solution  from 
one  str('amwis('  station  to  the  next.  For  more  details  of  the  practicalitic's  of  the  scIk'iik' 
th('  n'ader  is  reh’rn'd  to  Deni('r  k"  Hall  (1992). 

One  marked  difFer(’nc('  b('tw('<'n  our  pn'sent  work  and  that  of  D(’ni(T  k  Hall  (1992) 
is  that  w('  chose  to  calculate  the  nonliiwar  terms  in  physical  space  ratlu'i'  than  trans¬ 
form  sj)ac<'.  as  this  is  a  computationally  Vheaper'  policy  and  so  allowed  us  to  n'tain 
more  modes  in  our  calcuhitions.  This  necessitatc's  transforming  from  Fourier  spac(' 
to  physical  space,  ('ffecting  the  calculations  and  extracting  the  Fourier  coi'fficients.  It 
was  deci<'led  to  employ  Fast  Fourier  Tran.sforms  to  do  this  which  has  the’  beiK'fit  of 
n'ducing  the  cost  of  the  transformations  from  O(iV^).  (the  cost  of  n'ducing  A  nuxh's 
and  using  the  Fouri('r  transform  directly)  to  0(5.VlogA^  —  C.V).  TIk'  cod('  us('d  was 
based  on  the  original  Cooh-y  and  Tukc'y  (19G5)  algorithm  and  was  thus  limit('d  to  .V 
being  an  integral  power  of  2. 

Two  other  changers  wen'  made  to  tlu'  cod('  used  by  D('ni('r  «k:  Hall  (1992).  .4s 
was  d(’scrib('d  in  the  introduction,  tlx'se  authors  were  coiicernc'd  with  tlx'  (h'vt'lopnx'ut 
of  the  most  unstable  Giirth'r  mode  within  a  two  dimensional  boundary  layc'r.  To 
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achieve  this  most  unstalile  mode  Denier  fc  Hall  (1992)  took  an  arbitrary  distnr])ance 
and  then  marched  the  linear  stability  equations  for  a  long  distance  downstream  Ijefore 
changing  tci  the  nonlinear  system.  This  of  course  had  the  effect  of  ensuring  that  tin' 
most  unstable  component  of  the  original  disturbance  was  then  dominant  over  all  the 
other  components  so  that  when  the  full  nonlinear  equations  were  invoked  the  input 
vortex  flow  was  dominated  by  the  most  unstable  mode.  A  drawVjack  of  this  method 
is  that  the  lead-in  time  in  which  the  linear  mofle  develops  proves  to  be  a  significant 
l^roportion  of  the  total  computation  time.  To  alleviate  this  difficulty  we  used  the 
in'ogram  described  Ijy  Bassoni  i:  Otto  (1992)  to  compute  the  unstaltle  eigenfunction 
directly.  This  had  the  effect  of  making  the  lead-in  time  to  the  nonlinear  equations 
reduuflant  and  thus  we  could  start  our  computations  of  the  full  nonlinear  system  almost 
immediately.  Substantial  redtiction  in  the  computational  time  was  also  achieved  by 
using  a  stretched  grid  in  the  y  direction  (normal  to  the  cylinder  surface)  as  opposed 
to  a  \miform  one.  and  the  scheme  eventually  chosen  for  this  process  was  taken  from 
Macaraf'g.  Streett  ^  Hussaini  (19S8).  A  grid  is  requiretl  which  encompasses  the  region 
b('tw('en  (jp  =  0  and  ijp  =  Umax-  where  the  subscript  p  flenotes  the  physical  coordinate 
and  denotes  some  outer  bound  at  which  the  asymptotic  forms  of  the  solution  of 

the  system  (2.7)  are  supiKrsed  to  have  been  attained.  A  notional  computational  grid 
fl  S  Ur  <  1  is  introduced  which  is  related  to  t/p  via 

_  Umax  i^2  1)  l/r 

(<-2  ■ 

with  the  value  (cj  —  1)  taking  values  between  0.2  and  1.  The  quantity  (c2  —  1)  rep- 
rest'uts  tlu’  degree'  of  stre'tching  Ix'twee'n  the  large  and  small  steps  whilst  C]  controls 
the  rate'  of  stre'tchiug,  1  <  cj  <  G.  For  most  of  the  calculations  reported  here  we  used 
ci  =  2.4  anel  c^  =  2  which  allows  100  grid  points  to  be  distributed  between  Up  =  0 
and  Up  =  50.  The  re'sulting  distril>ution  e)f  peflnts  and  the  ce^rre'sponding  step  lengths 
are'  illustrate'd  in  figure  1.  Ne)tice  that  as  //  —y  t/mar  llit"  lengths  become  greater 
than  unity  which  in  a  finite  diffe’reuce  sche’ine  may  induce  errors.  In  the  pre.senit  case 
gradie'uts  of  the  elepe'neU'ut  functions  in  this  regime  are  minusc\ile  and  so  this  possible 
diffie'ifity  did  not  arise'.  This  particidar  non  luiiform  stretching  of  the  grid  points  al- 
lowe'fl  a  twe'iity  folel  reduction  iii  the  number  e)f  grid  points  over  the  regular  grid  used 
ill  Basseim  Hall  ( 1991 )  and  a  four  folel  improvement  over  the  piecewise  constant  steji 
le'ugth  grid  u.sed  in  Bassom  k:  Otte)  (1992). 

As  ])revie)usly  mentione’d  the  sefiution  strate'gy  u.se'd  lu're  is  essentially  that  de- 
scribe'd  in  Hall  (1988).  Suppose  that  the  solution  (F’n.  f))- c^n)  i^'^  known  at  some' 
spe'cifie'fl  statieiu.  .r.  and  supjieise  furthe'r  that  a  giu'ss  is  made'  for  the  solutieni  at  .r  -f-  e, 
^'(0)  The'  neinliue'ar  te-rms  on  the  right  hand  siele's  of  equations  (2.7  a.b) 
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were  calculated  using  these  guesses  and  the  system  solved  to  provide  updated  values 
These  updated  variations  were  then  put  into  the  nonlinear  t('rins 
and  the  process  repeated.  When  the  difference  between  successive  iterates  was  than 
some  norm  convergence  was  deemed  to  have  occurred  and  this  process  marclu'd  onto 
the  next  step.  For  a  computational  grid  composed  of  AI  points  and  a  calcidation  with 
N  modes  retained  then  the  convergence  norm  used  was 

*;=i\72-l 

E 

k=-.Nl2  3=1 

where  denotes  the  /;**'  iterate  of  the  A’**'  Fourier  comi^onent  of  the  flow  rpiantity 

U  evaluated  at  the  y-])osition. 

The  equations  (2.7  a,b)  may  be  discretised  for  any  particular  mode  into  the  forms 


|j-r(ll+l)  rT(«)| 

J.k  ^  \ 


A  iT/b'+i) 

I  ^  j.k  ^  j.k  I 


'^7(i  !’i71  +  2  +  ^7n*’r7i+l  +  <"771 1’777  +  <^777  ^’777  -  1  m  —2  +  .^777  ^  777  —  ^  777  - 

and 


(77?  =  -  1) 

(3.177) 


5^  777  *777+1  T  777  1/ 777  T  *  777  **  777  —  1  ■1"J  777  *'777  —  f  777  •  (  1)  (  3 . 1  (?  ) 


This  system  was  solved  using  a  technique  outlined  in  Appendix  B  of  Bassom  Otto 
(1992).  This  solves  the  probknii  (3.1)  in  one  sweep  ratlmr  than  treating  the  pair  (3.1 
a,b)  as  distinct  pentadiagonal  and  tridiagonal  problems.  In  essence,  the  equations  are 
solved  by  performing  an  outward  forward  elimination  followed  by  a  back  substitution. 

Earlier  we  alluded  to  the  fact  that  significant  computational  saving  was  obtained 
by  utilising  Fast  Fourier  Transforms  in  preference  to  decomposing  the  flow  quantities 
into  their  respective  components.  Standard  procedures  were  implemented  so  that  to 
change  between  physical  and  transform  spaces  we  write 


5’ 7 


-ikZj 


J=1 


1]- 


wh7're  Cj  =  27r(^  —  1)/.V  for  j  6  and  where  (]uantiti(’s  with  tikh's  are  in  tin' 

transform  space  and  those  without  in  the  physical  space.  Denier  ^  Hall  (1992)  wcuh' 
able  to  restrict  themselves  to  using  a  cosine  l>asi.s  due  to  the  nature  (if  Gdrtler  vortic('s 
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with  two  dimensional  l)oun(lary  layers  but  th('  addition  of  temporal  i)eriodieity  and 
crossflow  prevented  ns  from  doing  this. 

For  th('  majority  of  the  calcnlations  presented  htdow  the  nunilx'r  of  modes  rc'taiiu'd 
was  16  and  a  fairly  large  step  ns(xl  in  the  streainwise  flirection,  tj’picallv  f  =  0.01. 
Denier  &  Hall  (1992)  found  in  their  calcnlations  that  these  parameter  choices  gav(' 
results  to  within  graphical  accnracy.  Ftirther  testing,  using  100  jjoints  in  th('  normal 
(hn'ction  and  with  \jmai  =  50  has  confirmed  that  similar  choice’s  are  satisfactory  for 
the'  three’  elinie’nsie)nal  cemipntatienis  performeel  iiere. 


§4  Results 

We'  have’  eletailexl  the  nnmenie-al  methoel  by  whie’h  we  inve'stigate’el  the  se)lntie)n 
pre)i)e'rtie’s  eif  syste’in  (2.7)  althongh  a.s  yet  we  have  left  the  ele’finitiem  e:»f  the  amplitnele 
e)f  the  initial  vorte’x  nnspecifieel.  As  el<^scril)e<l,  we’  initiateel  onr  cennj^ntations  with  a 
nmltiple’  e)f  the  e’igenfnnctie)n.s  eef  the>  linearise'el  ve'rsie)n.s  e)f  system  (2.7).  Fe)r  a  spe'cifieel 
vorte-x  wave'immbe’r  k.  fre’epu’iicy  Q  anel  cre).s.sfle)w  A  the  methe)el  outlined  ley  Bassom 
k  Otte)  (1992)  was  nse’el  tee  ceempnte  the  e-eerre’sjeoneling  linearised  growth  rates  3^  anel 
the'  re'spective'  e'ige'ufnnctions  neermalise’d  .so  that  the  e’nengy  elefineel  ley 


is  eepial  tee  We'  re’fe'r  tee  A  as  the  amplitnele  eef  the’  initial  e'eendition. 

Onr  first  calcnlatieens  we're’  pe’rfeerme’el  jerimarily  as  a  ve’rification  eef  onr  ceeele  against 
the’  e’stablislie’el  re'snlts  eef  Denie’r  A:  Hall  (1992).  The'iie’e’  we’  cemsiele’reel  the  case  eef  a 
lenre'ly  twee-elinu’nsieenal  beennelary  laye’r  anel  steady  veertices  (A  =  =  0).  De’iiie’r 

Hall  (1992)  feennel  that  as  the’  veerte’X  e’veelve’s  ne)Tdine’arly  eleewnstre’am  tlie’ii  at  seeme’ 
peeint  the’ir  e’eemjentatieens  breeke’  eleewn  in  a  singnlarity.  Care’hd  inve’stigatieen  re’veale’el 
tliat  at  this  lex-atieen  the’  skin  frictieen  was  eef  the  fleew  vanishe’s  anel  the’ii  any  seelntion 
se’he’ine’  which  relie’s  njeon  a  mare’hing  teclmiepie’  be'cenne’s  invalielateel.  The’  vanishing  eef 
the’  skin  frie’tieni  was  inte’ipre’te’el  as  be’ing  inelicative’  e)f  the’  vejrtie-e’s  lere’aking  away  fre)m 
the’  wall  anel  me)ving  intee  the  e-eere’  eef  the’  benmelary  layer.  Onr  resnlts  for  various  iiiidal 
amplitnele’s  A  are’  inelie’ate’el  eni  figure’  2  whe’ie  we’  shenv  the  locatiem  of  the’  bre’akele)wn 
peeint  /  as  a  fnne’tie)n  e)f  A  fe)r  the’  niejst  nnstable’  GeArtle’r  me)ele’  with  ve)rtex  wave’iinmbe’r 
k  —  0.4  iO  (the  elistnrbane’e  was  introelne-e’el  at  .r  =  1  sei  the  distance  trave'lleel  ley  the’ 
pe’itnrbatieen  be’feere’  bre’akeleewn  is  ./•/,  —  1).  Xeit  surprisingly,  the  breakeleewn  leee’atieen 
.I'h  is  a  nuHieeteene’  eh'e’re’asing  fnnctieni  eef  the  amplitnele  A  anel  a.s  A  — >  0  see  Xf,  — >  cc 
•anel  the’  line’ar  preeble’in  is  re’trie’veel,  Feer  all  the  eether  e-ale'nlatieens  repeerte’el  npeen  he’re’ 
a  similar  re’iiel  is  eebse’rve’el  see  that  in  all  case’s  we  edieese’  A  =  0.2;  this  se’le’e’tieen  was 
made’  pnre'l}'  feer  illnstrati\-e’  pnrjeeeses  anel  has  nee  spe’cial  signifie’ance’  whatseeeve’r. 
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Ill  figurt'  3  wr  (leiuonstrati'  another  feature  of  the  two  -dinu'nsional  results.  Figure 
3a  indicatt's  the  dependeuee  of  the  lin  'ar  vortex  growth  ratt'  dr  hs  a  function  of  the 
vorti'x  wavenunilier  k  for  the'  first  two  modes.  As  is  now  well  docuinentc'd  hy  Denier 
rt  al.  (1991),  Denier  &z  Hall  (1992)  the  most  unstalile  modi'  has  growth  rate  dr  = 
0.312  at  k  =  0.476  and  dr  — >  0  as  A’  — >  0  or  k  — >  oo.  Wi'  note  howi'vi'r,  that  for  a 
significant  range  of  scaled  wavenuml)er.s,  roughly  0.3  <  k  <  O.S,  tlu'  vortex  growth 
rate  is  not  much  ri'duced  from  the  maximum  in  as  much  as  it  is  at  h'ast  00%  of 
the  maximum.  This  will  have'  imjiortant  consequences  for  later  results.  Figure  31) 
illustrates  breakdown  point  Xh  as  a  function  of  k  for  the  first  two  modes.  As  is  to  Ix' 
exjiected  for  fixed  initial  amplitudi'  A  and  fixed  wavenumber  the  first  modi'  always 
breaks  away  from  the  wall  before  the  second  one  and  this  trend  continues  for  higher 
modes.  However,  it  is  also  observed  that  for  the  primary  mode,  xi,  is  a  monotone 
decreasing  function  of  k.  Therefore  for  any  fixed  initial  amplitude  it  is  not  the  most 
unstable  mode  which  breaks  up  first.  Indeed  we  previously  remarked  that  for  a  whole 
range  of  wavenumbers  surrounding  the  most  unstable  value  linear  growth  rates  are  not 
too  different  from  the  largest  growth  rate.  This  has  potentially  im])ortant  consequences 
for  a  number  of  practical  flows  as  it  demonstrates  that  the  breakdown  proi)erties  of  the 
flow  are  crucially  dependent  on  the  nature  of  the  physical  characteristics  of  the  flow 
and  are  sensitive  to  the  nature  of  the  evolution  of  the  flow.  Alore  i)recisely,  suppose 
that  the  vortex  motion  starts  with  extremely  small  amplitude.  Th('n  one  would  expc'ct 
there  to  be  a  large  distance  over  which  the  motion  develops  essentially  in  a  linear 
manner.  During  this  time  the  most  unstable  linear  mode  would  overwhelm  modes  of 
other  wavelengths  and  once  the  vortices  had  grown  sufficiently  so  that  nonliiu'arity  is 
important  the  flow  l)ehaviour  would  be  dominated  by  that  of  the  most  unstable  mode. 
On  the  other  hand  if  the  initial  vortex  amplitufle  is  not  tiny,  nonlinear  effect  are  lik('ly 
to  b('  important  a  relatively  short  distance  downstream  l)y  which  stage  it  is  unlikely 
that  the  most  unstable  mode  dominates  the  others.  In  this  case  the  behaviour  of  tlu' 
most  unstable  vortex  is  not  likely  to  dictate  the  properties  of  the  breakdown  of  the 
flow. 

We  turn  now  to  consider  cp,ses  with  non-zero  crossflows.  As  noted  by  Bassom  L’ 
Hall  (1991)  we  may  restrict  our  attention  to  cases  in  which  the  crossflow  parameter 
A  >  0  since  by  suitably  transforming  the  systc'in  (2.7)  we  can  relate  flows  with  A  <  0 
to  appropriate  counterparts  with  A  >  0.  In  figure  4  we  recall  the  results  of  increasing 
crossflow  on  tlie  linear,  stationary  vortex  mode.  As  discussed  in  detail  by  Bassom  dc 
Hall  ( 1991)  the  effect  of  crossflow  on  linear  vortex  structures  of  wavelength  0(G~^'  )  is 
primarily  a  stabilising  on<'.  Indeed  figure  4  illustrates  that  when  A  >  0.410  the  vortc'x 
mode  is  stabilised  for  all  wavenumbers  in  the  (7(G^ )  regime. 
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Figures  5  a<5c  o  show  the-  variations  of  the  niaxiinuiii  growtli  rate  and 

the  eorr('sponding  vortex  wavenninber  k'mar  fnr  increasing  crossflow  parameter  A.  As 
already  noted,  increasing  A  lowers  and  it  is  the  mode  with  wavemiinher  h  ~ 

0.407  which  is  the  last  to  be  stabilised.  Further  k„iai  is  not  a  monotone  function  of  A 
as  might  have  been  anticipated.  Figure  5b,  as  well  as  showing  indicates  the 

growth  rate  of  the  inorle  with  wavenumber  k  =  0.470  (tin'  value  of  k,„ar  for  the  case  of 
zero  crossflow).  It  can  be  seen  that  tin'  difTerenc«“  in  growth  rates  of  the  most  unstable 
mode  and  that  with  k  =  0.47C  is  very  small  over  the  whole  range  of  A  considered. 
This  reinforces  the  earlier  comment  that  there  is  a  signiflrant  wavenumber  regime  over 
which  the  linear  growth  rate  of  the  vortex  nuxle  is  almost  constant. 

For  each  crossflow  A  we  integrated  the  aj^iiropriate  most  unstable  linear  eigenmode 
with  initial  amplitude  A  =  0.2  until  breakdown  occurred,  and  then  repeated  the 
experiment  with  vortex  wavenumber  k  =  0.476.  Th(“  results  are  summarised  in  figure 
5c.  It  is  seen  that  as  well  as  exercising  a  stabilising  influence  on  the  linear  mode, 
increasing  the  crossflow  tends  to  have  a  similar  effect  on  the  nonlinear  inodes.  Clearly 
the  greater  A,  the  greater  the  delay  downstream  before  the  initially  most  unstable 
mode  lireaks  away  from  the  wall.  Additionally,  for  the  same  initial  amplitude  A  and 
crossflow  A  the  mode  with  k  =  0.476  breaks  up  before  the  linearly  most  unstable  mode. 
Further  investigations  have  suggested  that  for  all  choices  of  A  and  vortex  frequency 
Q  then  if  two  linear  modes  of  different  modes  Aq  >  A’2  are  such  that  if  initially  have 
the  same  amplitude  then  if  they  are  marched  downstream  then  it  is  the  one  with  the 
higher  wavenumber  which  breaks  down  first. 

Finally,  we  rliscuss  in  more  detail  the  influence  of  unsteadiness  on  our  findings.  In 
tlu'ir  study  Bassom  &  Hall  (1991)  made  some  comments  concerning  the  properties  of 
tiiiK’-  dependent  linearised  vortices.  For  the  majority  of  their  work  these  authors  were 
jirimarily  interested  in  examining  neutrally  stable  modes  although  they  did  compute 
a  h'w  non-neutral  ones  (see  their  figvire  16).  Bassom  &:  Otto  (1992)  showed  that  for 
fixed  wavenumber  k  then  as  the  frequency  kl  of  the  mode  increases  .so  the  crossflow 
reejuired  to  maintain  neutral  stability  grows.  In  particular  it  was  shown  that  for  a 
non  dimensional  vortex  frequency  kl  the  stability  properties  of  the  vortex  are  sensitive 
to  th(’  sign  of  k}.  In  figures  6  a&  b  we  illustrate  a  facet  of  this  sensitivity.  For  each 
frecpK'iicy  il  and  cixissflow  parameter  A  we  show  the  wavenumber  of  the  most  unstable 
linear  mode.  k,„„j-,  together  with  the  growth  rate  of  that  mode  (/^r),,,^^..  For  <  0,  it 
it  is  obs('rved  that  as  A  increases  from  zero  so  A',„„j  decreases  whereas  for  >  0  this 
is  not  th('  case.  Correspondingly,  when  <  0  the  growth  rate  of  the  most  unstable 
mode  is  inonotoiK'  decreasing  in  A  whereas  when  S2  >  0  then  as  A  increases  from  zero  so 
there  is  a  crossflow  range  over  which  increa.ses.  This  increase  is  not  indefinite 

howev('r  ami  there  is  a  critical  A,  dependent  upon  kl,  after  which  the  growth  rate 
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decreases.  We  also  note  that  if  flj  >  Q2  >  0  then  inax^  flr  <  niax^  /ir  ^^^2:  A^ 

so  that  over  all  Cl  and  all  A  >  0  the  mode  with  the  greatest  growth  rate  is  stationary 
and  exists  in  a  two-dimensional  boundary  layer. 

The  breakdown  characteristics  of  unsteady  flows  are  described  by  figure  6c.  For 
each  frequency  Cl  and  A  we  marched  the  linearly  most  unstable  mode  of  amplitude 
A  =  0.2  from  j  =  1.  We  can  observe  the  somewhat  conflicting  roles  played  by 


crossflow  and  frequency. 

In  the  main,  for  a  prescribed  O  increasing  A  delays  breakdown  whereas  for  fixed 
A  and  increasing  Cl  this  phenomenon  is  enhanced.  Notice,  however,  one  important 
feature  which  runs  against  this  general  rule.  For  positive  Cl  then  a  small  to  moderate 
crossflow  actually  tends  to  promote  breakdown  although  larger  crossflows  do  reverse 
this  effect.  An  attempt  was  made  to  verify  this  trend  by  considering  larger  values  of 
Cl  than  those  illustrated  in  figure  6.  However  problems  were  encountered  as  Cl  grew 
and  these  difficulties  can  be  attributed  to  a  number  of  causes.  Following  on  from  work 
elucidated  in  Bassom  &  Hall  (1991)  it  is  the  case  that  for  small  crossflow  the  most 
unstable  linear  mode  first  has  a  small  wavenumber  relative  to  the  implied  scaling.  As  A 
grows  then  the  most  unstable  mode  corresponds  to  an  eigenfunction  that  moves  away 
from  the  wall  at  y  =  0.  At  Q  greater  than  about  2  this  movement  occurs  quickly  for 
small  changes  in  A  so  that  for  quite  moderate  values  of  A  the  eiegnmode  is  far  removed 
from  the  boundary.  As  found  both  by  Bassom  Hall  (1991)  and  Ba.s.som  Ss  Otto 
(1992)  the  numerical  solution  of  the  governing  equations  becomes  non-trivial  as  the 
vortex  moves  out  since  boundary  conditions  need  to  be  imposed  at  the  wall.  Clearly 
for  modes  concentrated  away  from  the  wall  large  changes  in  /i,  A  or  Q  can  lead  to 
almost  imperceptible  changes  in  the  values  of  the  eigenfunctions  at  the  wall  and  thus 
reliable  numerical  convergence  is  rendered  very  difficult.  However,  our  limited  further 
computations  for  Q  >  1.5  are  in  agreement  with  the  general  behaviour  described  above. 


§5  Conclusions  &:  Discussion 

In  this  work  we  have  detailed  the  nonlinear  spatial  evolution  of  unstable  Gortler 
modes  in  a  three-dimensional  boundary  layer.  In  particular,  the  roles  played  by  vortex 
wavenumber,  frequency,  and  the  crossflow  component  of  the  underlying  base  flow  have 
been  described.  We  feel  that  of  particular  importance  is  our  finding  that  (  all  other 
factors  being  equal)  of  two  modes  of  wavenumber  within  the  O(G^)  regime  the  one 
with  the  smaller  wavelength  will  be  the  first  to  breakdown.  This  then  suggests  that  in 
practical  situations  it  may  not  be  the  most  unstable  linear  mode  which  is  of  ultimate 
importance. 

In  many  cases  the  relevant  calculation  to  describe  the  breakdown  of  a  flow  is  oni' 
of  a  receptivity  type.  In  this  scenario  small  disturbances  within  the  boundary  layt'i 
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or  on  tli('  wall  of  tlu'  cylinclc'r  can  trigg('r  Gortlcr  vortices  and  tlu'  j)rccis('  method 
of  this  triggering;  fre(|neihly  ('xcites  mode's  of  a  preferre’d  wave’numhi'r.  If  this  occurs 
in  practice  then  our  calculations  ju'ovide  a  descrij)tion  of  the  evolution  of  tlu'  mode. 
Converst'ly,  if  a  range  of  wave'numher  modes  is  e'xcited  two  <'ventualiti('s  would  se'e'in  to 
he  j)ossil)l('.  First,  su])po.se  that  the  initial  distTirhance  is  vc'iy  small.  Then  it  is  to  Ix' 
e'xpe’ctt'd  that  the  pe'rturhation  trave'ls  a  long  way  downstn'am  before’  ne)nline’arity  has 
significant  e’ffect.  the  most  unstable'  line'ar  moele  will  be’  eleuninant  be’fe)re  this  peunt  anel 
its  e'vedutiem  characteristics  will  e’sse'iitially  ele’seribe  that  of  the’  whede  flenv.  Se’cemel. 
suj)])ei.se’  that  the  initial  perturbatie)n  is  ne)t  se)  small.  Our  results  summariseel  by  figure' 
3  have  shown  that  althe)ugh  the’ie’  is  a  unique’  me)st  tmstable  moele’  fe)r  e’aedi  e're)s.sfle)W  A 
anel  freepiency  Q  vertices  with  wave’uumbers  in  a  fairly  large  regieni  surremneling  that 
e)f  the’  me)st  unstable  moele  have  growth  rates  neit  ve’ry  elifferent  fre)in  the  maximum. 
The’refeue’,  l)y  the'  time’  nemlinearity  is  significant  it  is  imt  clear  that  the  nie)st  unstable’ 
nmele’  we)ulel  be’  ele)minant  anel  the’  bre’akele)wn  e’haracteristics  e)f  the’  whe)le’  fle)w  we)ulel 
invedve’  calculatie)ns  me)re’  inve)lveel  than  the)se’  repewte’el  he're’.  He)we've’r  we’  have’  she)wn 
that  fe)r  given  inerele’  amiditeiele  it  is  the’  liighe’r-wave’iiumbe'r  nie)eles  which  apj)e’ar  te) 
bre’akele)wn  first  se)  that  the’se'  ce)mpe)nents  of  a  spe’ctrurn  e)f  exciteel  me)eles  may  well 
pre)ve  te)  be  the  impe)rtant  e)ne’s. 

Denier  Sc  Hall  (1992)  she)we’d  that  whe’ii  the’ir  e'alculatie)ns  fe)r  ne)nliiiear  nie)ele's  in 
twe)-elimensie)nal  bounelary  layers  te’rminateel,  this  e'orrespe)nels  to  the'  skin  frictie)u  e)f 
the’  fle)w  vanishing  at  se)me  pe)int.  Once  this  happems  marching  scheme’s  as  useel  both 
lu’re’  anel  in  De’iiie’r  Sc  Hall  (1992)  canne)t  be’  continueel.  We  ceMifinn  this  fineling  fe)r 
e)ur  thre’e’  elime’nsie)nal  case’s  as  we’ll  but  we’  alse)  e)bserveel  that  beTore  the  skin  frictie)n 
vanishes  the  vele)city  pre)file’  ele'vele)p  inflectiern  pe)ints  at  pe)sitions  away  from  the  evall. 
The’  aj)i)earance’  e)f  the)se  inflection  i)e)ints  suggest  that  the  fle)W  will  become  susce'ptible 
te)  Rayle’igh  wave’s  whiedi  we)ulel  give  an  alternative  route  to  the  ultimate  breakeleiwn. 
The  analysis  e)f  these’  me)ele’s  woulel  be  of  inte’rest. 

Finally,  we  recall  that  all  e)ur  caloilations  have  be’e’ii  concerneel  with  consielering 
the’  evolutie)!!  of  j)erturbatie)ns  of  a  sjx’cified  wavenumber.  Of  course  in  se)me  situations 
a.  sj)e’e’trum  of  moeles  may  we’ll  be  present.  We  have  identifieel  situations  in  which  we 
might  expect  one’  mode  to  dominate  the  others  be’fore  nonlinearity  sets  in  but  in  the 
othe’r  case’s  calcidations  would  be  ne'e’fled  which  account  for  an  initial  perturbation 
which  contains  a  number  of  mode’s.  The  development  of  a  code  to  j)erforin  s\ich  calcu¬ 
lations  might  well  be  fe)rmidal)le  but  it  would  give  the  definitive  theoretical  description 
of  nonlinear  Geirtler  vortex  behaviour  in  thre’e  elime’nsional  boundary  layers. 
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Figure  1:  Distribution  of  grid  points  and  corresponding  sparing  for  the  ralndatious 
performed  here.  The  computational  coordinate  t/r  satisfies  0  <  y,-  <  1  an<l  is  r<'lat<'d 
to  the  physical  coordinate  ijp  by  i/p  =  b0xjr/{2  —  ■ 
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Figure  3: 

a)  Linear  vortex  growth  rates  lir  as  a  funrtion  of  wavenuinlx'r  k  for  the  two  most 
dangerous  modes  in  a  two  dimensional  boundary  lay<'r. 

b)  Breakdown  point  jt),  as  a  funrtion  of  wavenuml)er  k  for  the  modes  in  figure  3a  witli 
assumed  initial  amplitude  A  =  0.2.  The  symbols  flenote  rorresi)onding  points  on  the 
curves  of  a)  and  1)). 


Fip;iirf’  4:  Growth  rate  i  f  6r  of  the  most  dangerous  stationary  vortex  modes  within 
a  houiulary  layer  with  inrreasing  crossflow  A.  The  curves  correspond  to  A  varying 
between  0.005  and  0.410  in  9  equal  steps.  The  indicated  points  show  the  locus  of  the 
wavenumber  of  the  most  unstable  mode  as  the  crossflow  varies. 


n  )  \Vav('iimnl)<'r  of  the  most  unstable  stationary  mode  and  b)  its  rorr('s])ondin.c;  growtli 
I  'ltf'  as  a  function  of  crossflow  A. 

c)  Br('akdown  point  xh  as  a  function  of  crossflow  A  for  t)  tin'  most  unstabb'  mod('  and 
n)  the  mod<'  with  h  =  0.47G  for  initial  amplitmle  A  =  0.2. 
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